A Universal Criterion of Melting 
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Melting is analyzed dynamically as a problem of localization at a liquid-solid interface. A 
Lindemann-like criterion of melting is derived in terms of particular vibrational amplitudes, which 
turn out to equal a universal quotient (about one-tenth) of the molecular spacing, at the inter- 
face. The near universality of the Lindemann ratio apparently arises owing to strongly overdamped 
dynamics near melting, and despite the anharmonic interactions' being system-specific. A similar 
criterion is derived for structural displacements in the bulk of the solid, in particular the premelted 
layer; the criterion is no longer strictly universal, but still depends only on the harmonic properties 
of the solid. We further compute the dependence of the magnitude of the elemental molecular trans- 
lations, in deeply supercooled fluids, on the temperature and the high frequency elastic constants. 
We show explicitly that the surface tension between distinct liquid states, near the glass transition 
of a supercooled liquid, is nearly evenly split between entropic and energetic contributions. 



I. INTRODUCTION: WHAT IS "MELTING"? 

The word "melting" turns out to mean more than one 
thing. (For reviews, see dlEIESIMIESIES E3, E3, 

mo 

At a high enough temperature, the shear modulus of 
aperfect crystal would vanish, which, according to Born 
0, would lead to melting. This type of melting is 
sometimes called mechanical melting. Experience shows 
nevertheless that in real crystals, the lattice always be- 
gins to desintegrate at significantly lower temperatures, 
when the shear resistance of the bulk still exceeds at least 
a half of its zero temperature value P, 0, 13 • This type of 
melting, called thermodynamic melting, appears to usu- 
ally initiate at defects, most notably at the free surface, 
but also at grain boundaries, impurity sites etc. 
Another well known peculiarity of the melting/fusion 
transition is, while liquids arc relatively easy to over- 
cool, overheating a crystal seems very difficult, except 
under special circumstances, such as when the free sur- 
face is "clamped" , or in the case of water, which con- 
tracts upon melting thus allowing for internal melting 
Surface melting phenomena are very complicated, 
as is just about anything about surfaces: The melt may 
or may not wet the crystal surface; the thickness of the 
"premelted" layer, separating the crystal from the va- 
por, often exhibits power law scaling with the proximity 
to the melting temperature; the surface itself is usually 
reconstructed and often experiences roughening transi- 
tions below melting 0, [n|. One should add here the 
effects of lattice anisotropy 0, ^| and possible poly- 
morphic transitions near the melting temperature |16| . 
To summarize, conceptualizing melting as a dynamical 
process is not straightforward and depends on specific 
circumstances. 

In view of the complications above, it comes as a sur- 
prise that there should be the following simple, nearly 
univer sal p henomenological criterion, due to Lindemann 
[H El llSf : At melting the typical vibrational displace- 
ment, for a given crystalline class, should be some fixed 
fraction of the lattice spacing. Lindemann posited that 



the ratio should be about a half, implying direct colli- 
sions between the atoms constituting the lattice would 
become possible, leading to the lattice's demise. One 
may note that the Lindemann's argument accounts for 
anharmonicities in the problem but in a very generic fash- 
ion, through the existence of collisions. The value of the 
Lindemann ratio was later revised by Gilvarry [l9l | to be 
about one tenth and works rather well indeed: Data com- 
pilations 8, 20] show only a variation of 10% or so within 
a given crystal structure type, and the overall range is 
between 0.068 and 0.114. 

It is in the context of the Einstein's theory of vibra- 
tions in solids that Lindemann formulated his criterion; 
perhaps for this reason, the criterion has been a bench- 
mark in density functional studies (DFT) of crystal melt- 
ing (see e.g. j^H^), which are accurate in the short- 
wavelength, Einstein limit. Nevertheless, it was not un- 
til a DFT study of aperiodic crystals, when the Linde- 
mann ratio turned up on a first principles basis, as an 
order parameter: Wolynes and coworkers [23II24I ] demon- 
strated that a liquid, if failed to crystallize, should set- 
tle into (long- living) aperiodic structures. The transi- 
tion, or rather a cross-over, is characterized by a discon- 
tinuous change in the localization length from an effec- 
tively infinite to a finite value; the latter gives the vi- 
brational displacement at the mechanical stability edge. 
This length therefore directly corresponds to the Linde- 
mann length; its DFT-computed value matches well that 
observed in crystals, and, in the first place, the neutron 
scattering data in supercooled liquids and frozen glasses 
[25L l2l| . Lindemann-like criteria also naturally arise in 
treatments of energy landscape models of protein fold- 
ing and collapse [27L l28l . mean-field models of the struc- 
tural glass transition |29i . and in vortex lattice melting 
in superconductors 30]. Even in the absence of a first 
principles justification, the Lindemann rule is often used 
on purely empirical grounds, evidently owing to its re- 
markable circumstantial consistency, and its simplicity. 
Applications range from vortex lattices in rotating Bose 
condensates |31| to estimating the native state entropy 
of a protein [32|. Generalized Lindemann criteria have 
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bee n ap plied to defect-induced amorphization of a crys- 
tal |33| . or melting in one-component plasmas (see |34| 
and references therein). An inverse Lindemann criterion 
has been suggested for crystallization j3^. An increas- 
ingly useful application of the Lindemann criterion is in 
molecular dynamics simulations (see e.g. EH), where 
definitive observation of melting is usually beyond cur- 
rent computational technology. 

Perhaps, the most immediate objection to the Linde- 
mann's criterion is that it involves characteristics of only 
one of the two phases coexisting at melting; a proper 
criterion, presumably, should compare some property of 
both phases. For example, in the absence of extensive 
defects, comparing the bulk free energies would be ade- 
quate. In addition, as we now understand, a proper anal- 
ysis of dynamical melting should proceed with reference 
to processes at the liquid/solid interface. The present 
work implements these two notions in the following way: 

First, in Subsection III Al we consider the escape of a 
molecule from the solid/liquid interface into the liquid. 
Here, two length scales of molecular motions will arise, 
whose ratio to the molecular spacing is universal at melt- 
ing. One length scale is the size of the metastable mini- 
mum harboring a molecule which is about to change its 
location (for example, to exit into the liquid); the other 
is the extent of the transition state during the exit. The 
relation of the two scales to the vibrational amplitudes 
proper will depend on the detailed morphology of the re- 
gion in question; so for example, the surface roughness or 
the specific crystal face will affect surface melting. Under 
most circumstances, nevertheless, the derived criterion 
will simply amount to a simple Lindcmann-like criterion. 

Second, we ask in Subsection III Bl what would be the 
magnitude of the (weakly activated) molecular displace- 
ments in the bulk of the solid, in the presence of alter- 
native structural states. For glass- forming substances, 
these displacements are actually present in the bulk of 
the material, if it is supercooled. Otherwise, alternative 
structural states are present only in a "premelted" layer, 
if any. A proper melting criterion will be formulated, 
which amounts to comparing the vibrational molecular 
amplitudes to the structural displacements: in a stable 
solid, the former should less than the latter. We will fur- 
ther deduce the dependence of the structural displace- 
ments on the material's stiffness, and the temperature. 
The results apply directly to premelted layers and super- 
cooled liquids. 



II. DERIVATION OF A MELTING CRITERION 

A. Melting at the Interface 

It will be most straightforward to see how a melting 
criterion arises for molecules that are directly at the solid- 
liquid interface, in the sense that here, simply an isolated 
activated event is required in order for a molecule to exit 
irreversibly into the liquid. A number of conventions as 



to what an "interface" is are possible and are subject 
to the same ambiguities as the definitions of the phases 
themselves. It is beyond reasonable doubt that the sur- 
face region of a melting solid is far from the simplified 
text-book pictures, even ignoring surface roughening ef- 
fects: The interface is not sharp, and often extends for up 
to several tens of atoms, as could be deduced, somewhat 
indirectly, from studies of (pre)melted layers demarcat- 
ing a solid from its vapor [Til 13^ . (See also a recent 
review [l2| on surface ice melting.) That the interface 
is "diffuse" near melting has been also concluded the- 
oretically, via density functional studies [3^ employing 
specific ansatzes for the equilibrium distribution func- 
tions in the solid and liquid. Furthermore, the hetero- 
geneity across the layer is not only structural but must 
also involve a heterogeneity in relaxation times, i.e. the 
lifetimes of long-living local structures. Similar to the 
density, the life-times in the premelt layer interpolate be- 
tween those in the liquid and the solid. Several ways to 
deal with the ambiguity in defining a solid-liquid inter- 
face may be proposed. For example, Oxtoby and Haymet 
|38| employ an appropriate equilibrium order parameter 
changing continuously when going from liquid to solid. 
Trayanov and Tosatti [39j analyze premelted layers, in a 
mean-filed fashion, based on two order parameters, "den- 
sity" and "crystallinity" . Here, in order to focus on the 
dynamical aspects of melting, we will use the life-times 
of local metastable structures to establish an operational 
criterion of whether a molecule is in the solid, or part of 
the liquid. While the equlibrium interfacial region may 
be discussed only in a broad sense, as a "diffuse" entity, 
the dynamic interface will turn out to be thin and well 
defined. 

We will distinguish between the liquid and the solid 
in the usual way, via symmetry, and will specifically fo- 
cus on the time scale on which the symmetry is bro- 
ken/restored. Consider a substance consisting of a sin- 
gle, relatively compact molecular species, in the classical 
regime. The crystal breaks the translational symmetry 
in that here, one can label the molecules based solely on 
their each being located within a particular, well defined 
cell. (C.f. however the incommensurate quantum crys- 
tals j43|). One may speak of a fluid, on the other hand, 
when such labelling is impossible. The corresponding 
translational symmetry is physically maintained by par- 
ticle transport. Call To the time it takes a molecule to 
diffuse a distance defining the volumetric density of the 
liquid. Choose a compact, specific cell, in space, whose 
volume is equal to the volume per molecule in the liq- 
uid. Since the time To is significantly longer than the 
time scale of density fluctuations, it is guaranteed that 
another, identical molecule will have visited the chosen 
cell, within time To upon the exit of the previous cell's 
dweller, thus erasing the possibility that one be able to 
label a molecule by its spatial location. It is therefore at 
times exceeding t that one may speak of a liquid state. 
Recall also that we are considering compact molecules, 
and so rotational diffusion/equilibration is not an issue. 
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It is useful, for future reference, to compute the particle 
exchange time t , in terms of the collisional time: Let 
1/a 3 = n be the molecular concentration in the fluid, so 
that a is the average, volumetric molecular spacing. The 
typical collisional, or auto-correlation time (also defining 
the density fluctuation time-scale) is T au to = m/C,, where 
Q ~ 6nar] is the friction coefficient, ry is the viscosity, and 
m is the molecule's mass. The time To it takes to diffuse 
a distance a is roughly a 2 /6D, where D is the diffusion 
constant, related to the friction through the Einstein's 
relation D = k B T/(. As a result, 

to _ 67rVa 

Tauto pksT 

where p is the liquid's mass density. One may directly 
check that near melting, this ratio is generically about 
10 3 but varies within an order of magnitude or so between 
different substances: For instance, cobalt and sodium 
yield 1.3 • 10 3 and 2.1 • 10 3 respectively. This indicates 
it takes about a thousand molecular collisions or so, per 
molecular volume, to establish local thermal equilibrium 
in a liquid. Finally, note that the large value of the ratio 
in Eq.Q is an internal test of the argument's consis- 
tency: To give a counter example, it would be incorrect 
to use analogous logic to estimate equilibration times in 
dilute gases. Here one would find, using the elementary 
kinetic theory, that To/T auto — 9(er/a 2 ) 2 < 1, where a is 
the molecular scattering cross section; clearly To does not 
correspond to an equilibration time scale. Of course, the 
rate limiting step during equilibration in dilute gases is 
diffusion in the momentum space, which is responsible 
for establishing the Maxwell distribution of velocities; 
whereas in dense liquids near fusion, the rate limiting 
step is configurational equilibration. 

Consider now a region of space occupied by a solid 
and its melt, at some temperature T just above the low- 
est temperature, T m , at which surface melting is possible. 
Suppose there is a molecule, in the region, that fails to 
move a distance a in the time To . (To avoid confusion we 
note that if long wave-length sound is present in the sys- 
tem, one should stipulate that the local reference frame 
move with the sound.) All such molecules can not be 
regarded as part of the liquid, as just discussed, and so 
we must regard them as part of the solid. The boundary 
of any (spatially) closed set of such molecules may there- 
fore be defined as the solid-liquid interface. On the other 
hand, the inability of a molecule to move the distance a 
in time To is equivalent to saying the molecule is residing 
in a metastable free energy minimum. In other words, 
a molecule is part of the solid, if the escape time T esca pe 
from its current neighborhood exceeds the exchange time 
To: 

^escape > t . (2) 

The escape time T cscapc will generally differ for distinct 
crystalline faces, or various distinct surface morphologies. 
The above stipulation that 

T = T m + 0+, (3) 



V(x) 




FIG. 1: A generic schematic of an escape free energy profile is 
shown. Its, the transition state size, demarcates the vicinity 
of the saddle point within the thermal energy from the top. 
d and dxs are defined in text. 

implies that there is at least one specific face/morphology 
which is in near equilibrium with the liquid, and there 
are no face/morphologies which arc melting in a sponta- 
neous fashion. In the following we will specifically con- 
sider those faces that are melting in a quasi-equilibrium 
fashion. For these, 

Tescape = To. (4) 

One may regard this expression as the dynamical def- 
inition of a solid/liquid interface. A melting tempera- 
ture T m , as defined above, will generally differ from the 
usual calorimetric melting temperature, if the solid melts 
anisotropically. (The latter is usually the case.) 

In order to estimate T eS cape> we shall adopt the ap- 
proach of Frauenfelder and Wolynes (FW) j^jj (see also 
[42.l43j l. who have delineated the various activated trans- 
port regimes, and computed the corresponding rates, in 
terms of several characteristic length scales. First note 
that owing to the frequent collisions (see also below), 
the motion of the reaction coordinate corresponding to 
the escape mode, is strongly overdamped. An opera- 
tional criterion of this is that the particle's mean free 
path be significantly shorter than the transition state 
size: i m f p -C Zts- (A schematic of the free energy profile 
along the progress coordinate is shown in Figd) As a re- 
sult, the particle stays in the transition state region for a 
long time, relative to the molecular collisional time T auto , 
leading to a large number of barrier crossings, while at 
the top of the barrier: N c = Zxs/^mfp- The corresponding 
rate may therefore be estimated using the standard tran- 
sition state result, multiplied by a (small) transmition 
factor k ~ 2Z m f p /Z T s = 2T auto wth/'TS, as appro priate in 
the overdamped Kramers limit. Here, Vth — y/BkbT/m 
is the thermal velocity of the particle. As a result, the 
Kramers rate reads: 

-X = 1 (H) gTauto^th -yt/k B T 
eseape 3 ^ ' 
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where 



d = 



/2tt 



dx exp 



-[V(x)-V(x m )]/k B T 



(6) 



is the "size" of the metastable basin. (The label x m in 
the integration limits indicates that the integral is un- 
derstood in its asymptotic sense in terms of the basin's 
width.) The coefficient l/y/2n was incorporated in 

1/2 

the definition so that d — ((Ax) 2 ) , if the poten- 
tial is strictly harmonic at the minimum x m . (\v\) — 
y^2/(37r)?;th is the thermally averaged particle speed that 
enters the expression for the molecular flux at the barrier 
top. Lastly, V* = V(x^) — V(x m ) is the barrier height. 
The latter is actually quite easy to evaluate after one real- 
izes that an escape event will have occured, if a particle's 
displacement just exceeds the typical thermal vibrational 
amplitude of the interface. This is because the particle 
will have crossed the thermally definable boundary of the 
solid and must be regarded as part of the liquid. For this 
simple argument to be valid, it is essential, again, that 
the specific face of the solid be wetted by the liquid, i.e. 
the face be actually melting. Since the typical energy of 
the surface vibrations is exactly ksT, the particle's typ- 
ical displacement free energy cost will not exceed ksT, 
hence V* = ksT. Further, we define 



d 



1 



TS 



2%/2 



(7) 



where the numerical factor is chosen so that if the bar- 
rier is parabolic at the top, then dxs = {{Ax) 2 ) in 
the inverted potential at the saddle point. By definition, 
to(w*) 2 (Zts/2) 2 /2 = ksT, where u* is the under-barrier 
frequency, see FigQ] After putting all this together 



and recalling that tq = a 2 /6D 
a 2 /2v 2 h T auto , one has by Eq.Q: 

ddTS 1 



a m/6fcsTr a 



4\/6 



0.01, 



tic: 



(8) 



universally. We therefore observe that it is possible to for- 
mulate a purely kinematic criterion of melting, in terms 
of the ratio of length scales characterizing molecular mo- 
tions in the interface region. The numerical constants 
on the r.h.s. of Eq.(jHJ should not, perhaps, be taken 
too seriously; nevertheless the estimate is expected to be 
accurate within a factor of 2 or so. 

A number of comments are due here. First of all, can 
one even apply a transition state argument, when the bar- 
rier is so low? The answer is yes, because of the high fric- 
tion. Indeed, suppose for a moment there were no barrier. 
Even so, the exit would be far from instantaneous, being 
subject to (frequent) collisions with the nearby molecules, 
just as are the molecular motions in the neighboring liq- 
uid. Using the Smoluchowski diffusion limited reaction 
rate expression (in ID), one still gets the same basic scal- 
ing for the escape rate: k cscapc ~ D/l^ s ~ -y 2 h r aut o/^TS 
(see also [12 )• The overdamped character of the molecu- 
lar motion is essential in the present context, and so one 



should like to estimate the actual value of the Z m f p /ZTS 
ratio, which is effectively the small parameter of theory. 
The estimate in Eq.Q suggests that this ratio is indeed 
quite small. Numerically it is of the order one hundredth, 
considering that dxs /a ^ 0.1. At any rate, molecular mo- 
tions in liquids are overdamped, near T m , and so must 
be the motions in the corresponding crystals, near T m , 
because they are of comparable (usually greater) density. 

Despite simplifications due to the high friction, the 
shallowness of the metastable potential complicates the 
interpretation of the simple result in Eq.JHJ). Since the 
transition state size is virtually equal to the reaction path 
length, it probably makes little sense to differentiate be- 
tween dxs and d. For those same reasons, the integral 
in Eq.© is only meaningful in the lowest order in the 
reactant basin width, the expansion being asymptotic of 
course. This leads to a simpler yet relation 



dL 
a 



0.1, 



(9) 



where dt ~ d ~ dxs stands for the amplitude of the 
reversible motion in the molecular metastable minimum 
at the liquid-solid interface. The subscript "L" alludes 
to "Lindemann" ; the parameter dh will be sometimes re- 
ferred to as the "Lindemann length" . Importantly, di is a 
measure of the displacement in the direction of the fastest 
escape, which is perpendicular to the surface. Further, 
a signifies the molecular spacing right at the interface. 
Strictly speaking this implies our (volumetric) molecular 
spacing is a variable changing continuously across the in- 
terface. This appears reasonable as the crystal and liquid 
density differ actually relatively little, not by much more 
than the typical thermal density fluctuations in the crys- 
tal, near the melting temperature. (See also the earlier 
mention of the Oxtoby and Haymet's order parameter 
|38f.1 Finally, the parameters d and dxs should be der- 
erminable in a simulation, with a modest computational 
effort. 

Are the simple relations JSJ and © consistent with 
the general notions of surface melting? Yes, in a rather 
plain way. To give an example, suppose the surface is 
rough, so that it has a corner, or an edge. Clearly the 
vibrational amplitudes at corners and edges are larger 
than those at extended flat faces, because there are fewer 
neighbors. This is consistent with the expectation that 
corners (and edges) melt first, i.e. at lower temperatures 
than say flat faces. The relation in Eq.J5J) is quantita- 
tively consistent with Valenta's calculations of the vibra- 
tional displacements at the three distinct crystall faces 
of lead, i.e. (110), (100), and (111) Q. These faces 
are known to premelt at increasingly higher tempera- 
tures. According to Ref.0], the vibrational displace- 
ments at all three faces turn out ot be of nearly the 
same magnitude of 0.22A or so, at the respective pre- 
mclting temperatures. In any event, one should expect 
that denser packed, stiffer faces will exhibit lesser vibra- 
tional displacements normal to the interface and there- 
fore will melt at increasingly higher temperatures. This 
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qualitative notion is consistent with the results of MD 
simulations on an FCC crystal, by Kojima and Susa [45| . 
and on a BCC crystal, by Sorkin at el. 0|. 



B. A thermodynamic melting criterion in the bulk 

Let us now discuss the implications of the relations (JSJ 
and jnj for the molecular motion amplitudes in the bulk of 
the solid. The typical vibrational displacements, d V ib(0), 
at the surface and in the bulk, d v ih{z), will be certainly 
comparable. The variable z measures the distance from 
the interface, see Fig'0 In the bulk, the greater lattice 
stiffness in the bulk will be partially compensated by a 
smaller lattice constant, save the substances that expand 
upon freezing. As a result, the local (d v ib(z) / a(z)) ratio, 
while also generically of about 0.1 in value in the bulk, 
is no longer expected to be strictly universal, in con- 
trast to that at the interface (see below) . Now, one must 
bear in mind that motions resulting in a locally different 
structure, may also be present on the solid side of the in- 
terface. In the preceding Subsection, we have computed 
the magnitude of the displacements dh such that would 
lead to the exit of an atom from the solid into the liq- 
uid. In this Section, we will compute the magnitude of 
analogous irreversible atomic displacements, but inside 
the solid. Here, the atom also exits its present lattice 
site but to find itself in a reconfigured lattice, not in the 
liquid. Given the lattice has reconfigured and is poised 
to accept a particle, the latter will transfer in a nearly 
activationless fashion, similarly to exiting into the liquid. 
The lattice reconfiguration itself is cooperative event con- 
sisting of a large number of those elemental, nearly acti- 
vationless transitions occuring on the time scale r eS cape- 
The latter must occur in a concerted fashion, implying 
the cooperative reconfigurations are rare and occur on 
much longer times scales: 



Tstr ~> r e: 



(10) 



where the index "str" indicates "structural". 

First of all, do such rare structural reconfigurations 
take place in the bulk? We argue in the following that 
they do indeed, within premelted layers. According to 
the surface calorimetric studies of Santucci at el. |46j . 
the excess entropy of the "premelted" surface layer of the 
Li (110) face, relative to the bulk crystal, is about a half 
of the bulk liquid entropy excess, see Fig. 3(b) of Ref.[46|. 
(The author is aware of surface calorimetric estimates 
only on this particular substance, however comparable 
values of the surface entropy are expected for other mate- 
rials as well, see and references therein.) This clearly 
implies that the excess entropy, s c , within the premelt 
layer, is intermediate between that of the crystal and the 
liquid. Furthermore, this entropy must decrease into the 
bulk, so as to interpolate between the liquid and the solid 
values. The latter is zero. The relatively high density of 
states in the premelt, corresponding to s c ~ fcs/2 per 
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FIG. 2: A cartoon of a solid-liquid interface is shown, in- 
cluding the premelt layer; the latter is characterized by some 
disorder. (Here, we ignore possible roughness or reconstruc- 
tion of the surface.) The area containing the doubled circles 
illustrates a region undergoing a cooperative structural recon- 
figuration. The cooperativity is necessary to minimize density 
variations during the transition. The solid and the dashed 
circles denote the atomic configurations before and after a 
structural transition. An individual, elemental displacement 
is of magnitude di(z) mildly dependent on the distance to the 
interface. Each of these elemental displacements are nearly 
activationless; whereas the total cooperative event has a high 
barrier because the probability of a large number of concerted 
elemental events is low. 



particle, would be impossible to account for by trans- 
lations of vacancies: The vacancy concentration would 
be too small, considering that the density and the com- 
pressibility of the premelt are comparable to those in the 
solid bulk. (The vacancy formation energies are in the 
eV range, of course, see e.g. 47].) As a result, the room 
for molecular translations is provided not by diffusing 
vacancies, but some other structural degrees of freedom 
that involve more than one particle. Another possibility 
is dislocations, whose significance seems less straightfor- 
ward to estimate than that of vacancies. At the present 
stage of theory, however, it appers that dislocations be- 
come important closer to the point of mechanical melting, 
which they most likely orchestrate in the first place. The 
mechanical melting temperature seems to be several tens 
of degrees above the bulk thermodynamic melting point 
(see [EH) 111 and references therein.) 

What is the detailed microscopic nature of the con- 
figurational degrees of freedom, in the surface premelted 
layer? It is likely that they correspond to transitions 
between long-living structures analogous to the local 
metastable structures predicted by the random first or- 
der transition (RFOT) theory of supercooled liquids and 
glasses, of Wolynes and coworkers 0>^ii|5(J- The tran- 
sitions span regions 3 to 6 molecular units across, and 
have been directly observed by a variety of non-linear 
spectroscopies, see [5lj | for a recent review. The config- 
urational entropy of a supercooled liquid, just above the 
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glass transition, is about 0.8ks per molecular unit. It is 
therefore comparable to that in the premelted layers. To 
complete the analogy, note the viscosity of supercooled 
melts, near T g , is extremely high, consistent with the pre- 
melted layers exhibiting progressively slower relaxation 
away from the interface proper. 

One may estimate, semi-quantitatively, how the mag- 
nitude of the elemental displacements decreases away 
from the surface, as the stiffness increases. The mag- 
nitude of the displacements will generally depend on the 
distance from the interface, i.e. d L = d^z), see Fig|21 
Because a particle exits not into the liquid but into a 
matrix poised to receive the particle, we no longer have 
the convenience of the barrier being equal to the thermal 
energy. We may say, nevertheless, that the barrier will in- 
crease into the bulk because the lattice stiffness increases 
away from the surface. Let us denote this z-dependent 
barrier as V$(z), so that 



V*(0) = k B T m . 
Now, it turns out that both T au to and r, 



(11) 



escape 



are z- 



independent, because of the detailed balance: There is 
no net particle transport at these time-scales. The ac- 
tual net particle transport occurs at the much slower time 
scale, r s tr, of the extended structural transitions, and is 
consistent with the existence of a density gradient in the 
premelt layer. (The detailed balance argument above is 
only valid when the elemental and the structural tran- 
sitions are time scale separated, which is indeed true, 
in view of Eq.^nj)). The constancy of r au to and r eS cape, 
together with Eqs.JSJ) and yield that the following 
quantity is invariant throughout the solid: 



d L (0) 



■ YHz) 

e v, w = e. 



(12) 



This statement can be used to self-consistently determine 
the value of the dh ( z )/dL(0) ratio, after one recalls that 
the activation barrier V* arises from the elastic strain 
of the lattice: Since the thermally relevant vibronic dis- 
placements are within about one-tenth of the molecular 
spacing or less, the elastic energy, as a function of the 
lattice strain, is dominated by the quadratic component, 
i.e.: 



V x {z) ~ K{z) 



d L (z) 



a(z) 



a\z) = K(z)a(z)dl(z), (13) 



Here, K(z) is a local effective elastic modulus that de- 
pends both on the isotropic compressibility and the shear 
modulus, and a(z) is the local volumetric spacing; in 
other words V*(z) gives the elastic energy arising from 
the local strain [d L (z)/a(z)] 2 , per molecular volume. 
(The ratio [c?^(z)/a(z)] 2 is a measure of the thermally 
averaged square of the dimcnsionlcss elastic strain used 
here as in the standard elasticity theory |Hj.) One ob- 
tains, as a result: 



dh{z) 



d L (0) 



(14) 




FIG. 3: Shown are the dependences of the Lindemann and the 
vibronic displacements, relative to their maximal value dmax, 
achieved at the surface, as functions of the dimensionless stiff- 
ness The latter is the product of the elastic constant and 
the lattice constant, in terms of the minimum value of this 
product achieved at the surface. 



where, again, [di(0)] ~ 0.1a and the parameter £(z) gives 
the stiffness of the lattice, at distance z from the inter- 
face, relative to its value at the interface: 



CO*) 



K(z)a(z) 
K(0)a(0) ' 



(15) 



Note that K(0)a(0)d 2 L (0) = V i {0) = k B T m . We stress 
again that the (very low) barriers, implied in Eas. (jl^j) 
and HI 4(1 , correspond to elemental translations that would 
occur given another structural state exists and therefore 
apply to all of the bulk. Because the elemental events 
are subject to the existence of an underlying structural 
transition in the corresponding region, it is appropriate 
to term the displacements g?l as fiduciary displacements, 
or presumable displacements. 

It is easy to solve, numerically, Ea. l(14|) for the 
dz,(z)/<iL(0) ratio as a function of the dimensionless stiff- 
ness C( z )j with the result shown in Fig|3J The fidu- 
ciary displacement d^ (z) should be compared to the ther- 
mal vibrational amplitude d v ib(z), which is fixed by the 
equipartition theorem: 



K(z)a(z)d 2 vih {z 



k R T. 



(16) 



The l.h.s. of the first equation gives the usual elastic en- 
ergy density of the solid, times the volume of a unit cell. 
As follows from the discussion in the preceding Subsec- 
tion, 



dvib(0) = d L (0) = d n 



(17) 



The vibronic displacement d v ib( z ) depends on the stiff- 
ness in a very simple way, and is shown in FigO alongside 
the C-dependence of the fiduciary, Lindemann diplace- 
ment. Note that here we did not have to treat explic- 
itly the anisotropy of elastic properties, as long as we 
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were consistent in measuring the displacements c?i(0) and 
^vib(O) normal to the free surface. 

Since the z-variable is a dummy label, it is appropriate 
to simply think of the displacements and d v ib ; rela- 
tive to their maximal values achieved at the surface, as 
functions of the dimensionless stiffness (. According to 
Fig|3J it is always true that 

d vih (z) < d L (z), (z > 0), (18) 

c.f. Eq.©. We therefore observe that solids may be 
defined as collections of molecules in which the magni- 
tude of purely vibrational excitations is always smaller 
than the extent of the fiduciary structural displacements. 
As a result, a solid is (typically) capable of sustain- 
ing vibrational excitations without irreversible structural 
changes. According to Fig|3 a stabler solid corresponds 
to greater values of the stiffness parameter £, correspond- 
ing to larger values of the difference (<iz, — d V ib)- One 
may therefore interpret Ea. H18|) as saying that the stiff- 
ness of a lattice that allows for alternative structures, is 
somewhat smaller than the stiffness of the correspond- 
ing lattice which is strictly elastic, i.e. mechanically sta- 
ble. We thus observe that vanishing of the shear modulus 
is not required in order for the lattice to be locally me- 
chanically unstable. Now, the limiting situation where 
dvib(z) — d>L{z), at z = 0, bears a dual meaning: (a) 
it signifies melting and (b) it defines the boundary be- 
tween the solid and the liquid phase. In view of Eg. (1121) . 
the statements in this paragraph are valid regardless of 
the elastic free energy being strictly quadratic in the dis- 
placement. Finally note that Eas. ll4fl and (|16|1 are laws 
of corresponding states for the structural and vibrational 
displacements in a solid. 

We have thus found that formulating a criterion of 
melting and defining a solid, in the first place, requires 
introducing fiduciary structural modes, in the sense that 
these modes may or may not be readily observed in a real 
material. We have argued that such modes probably exist 
in a sufficiently extended premelt layer. Otherwise, there 
is still a way to force the solid to sample the spectrum of 
those structural excitations, even at temperatures below 
the melting point, namely by quenching the correspond- 
ing substance from its liquid state into a frozen glass. 
Many polymers simply do not form crystals, and are dif- 
ficult to characterize morphologically in a definite way, 
in the first place. In these, there is no shortage of struc- 
tural transitions. A similar comment applies to folded 
proteins. Certain substances, that are very difficult to 
supercool, such as ice, amorphize under pressure (see e.g. 
[53j). Finally, energetic particle beams may be employed 
to locally force a crystal into a structurally excited state. 
The possibility of such externally induced structurally ex- 
cited states is actually a physical realization of a Maxwell 
construction! To give an analogy, the Maxwell construc- 
tion in a non-ideal gas also uses presumable states, i.e. 
such that (dp/dV)r > 0, in deriving a criterion of a first 
order transition, i.e. boiling/condensation. 

Now, suppose there is a significant premelt layer on 



the surface. Direct spectroscopic observations of these 
are intrinsically difficult, since a beam (neutron or X- 
ray) sufficiently intense to produce useful signal, would 
also inevitably heat up and liquefy the surface. On the 
other hand, surface, near-atomic resolution techniques, 
such as those employed by Israeloff [HiJ, may be helpful 
in characterizing structural surface modes. Importantly, 
a conclusive surface study must be non-linear so as to 
sense dynamic heterogeneity, in order to distinguish ir- 
reversible structural relaxations from purely vibrational 
excitations. 

Finally, for completeness, one should like to see 
whether the arguments above would also be robust in 
systems other than single-component liquids. Unfor- 
tunately, mixtures necessitate considering many addi- 
tional factors, well beyond the scope of this article, 
primarily owing to the details of interaction between 
the constituents of the mixture, phase separation etc. 
Here, we will limit ourselves to a few remarks. Suppose 
the molecules are (chemically identical) rigid, relatively 
weakly interacting rods. Clearly the melting will occur 
from the face parallel to the orientation of the rods, oth- 
erwise the molecule would have to slide out of the ma- 
trix, which would be too costly energetically. (We as- 
sume the rods are closely packed, in a colinear fashion.) 
The "melting displacement" d is therefore perpendicular 
to the rods' orientation, and, consequently, is also about 
one tenth of the monomer size, but probably not univer- 
sally. Analogous logic applies to a crystal made of weakly 
interacting disk-shaped molecules, which will melt off the 
crystallographic plane parallel to the disk planes. (This 
is in the case when the discs are stacked in a coplanar 
manner.) The case of mono-valent ionic melts is proba- 
bly adequately viewed as that of a single-component sub- 
stance with a slightly elongated molecule, since the melt 
will tend to be locally electrically neutral, on average. 
Again, we recover the basic estimate from Eq.© where 
a is the volumetric spacing between the distinct moieties, 
not the unit cell size; d is the displacement between the 
neighboring atoms. 



III. ACTIVATED MOTIONS IN SUPERCOOLED 
LIQUIDS AS LOCAL MELTING 

The analysis of Subsection III Bl was conducted on the 
assumption that solids exist, of course. In other words, 
the argument by itself could not be used to establish the 
existence of a solid state, but could only give the condi- 
tions of stability of a solid once it is formed. The details 
of solidification, below T m , are known to depend on how 
fast one cools the substance and on the liquid's viscos- 
ity. If either of the cooling speed or the viscosity is high 
enough, the nucleation of the crystalline phase will be 
prevented, with the liquid finding itself in a supercooled 
state for a significant amount of time. 

Here, we will estimate the temperature dependence of 
the elemental displacements in a liquid that was cooled 
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below its fusion point, but has failed to crystallize. Since 
the structural equilibration has never occured, displace- 
ments must take place such that they imply a change 
of the local structural state. Since these displacements 
are at the mechanical stability edge of the material, we 
may also denote their magnitude with c2l(T), where we 
explicitly indicate the temperature dependence of the 
magnitude of the elemental displacements. (Clearly, 
dL{T m ) = cLl{z = 0).) Therefore in a supercooled liquid, 
the displacements d^ are no longer fiduciary, but strictly 
factual. On the other hand, the regular liquid state, that 
does not discriminate between molecular displacements 
of less than size a (recall our "labelling" discussion) , now 
becomes a fiduciary state! (For reference, we point out 
such a regular liquid state explicitly arises, for instance, 
in the DFT study of aperiodic crystals in Ref. 23], or of 
the mean- field Potts glasses |5q.) We may use this new 
fiduciary state to write the law of corresponding states 
from Ea. (|14|l at a temperature T below T m : 



d L (z,T) 



1 2 



d L (0,T) 



(19) 



Here, even though the z-variable is dummy, it is strictly 
implied that z — oo, i.e. the interface is infinitely remote 
from the bulk, by construction. Since the universal quan- 
tity dL = dr,(0)/a(0) is, by its very meaning, tempera- 
ture independent, it is convenient to rewrite the above 
relation in the following form: 



(d L /a ) 



2 i{Ka 3 /T) <AiJ^_ 

e l 



where 



((Ka'/T) ee 



K(T)a 3 (T)/k B T 



K(T )a 3 (T )/k B To' 



(20) 



(21) 



As a reminder, K (T )a 3 (T )d 2 / k B T = 1, for any T . 
While, in principle, the interface at any temperature 
Tq < T m may be used as a reference state, in practice 
one can measure the stiffness only at the physical melting 
temperature T m ; therefore most conveniently one would 
set To = T m in Ea. H2l)|) . Again, we have arrived at a 
law of corresponding states, even though the variables 
are distinct from those entering Ea. H14fl : The displace- 
ments are in relation to the lattice spacing; the running 
coupling constant Ka /T, governing the displacements 
magnitude, is dimensionless too, and is of the order 1. 
The latter suggests that below the melting temperature 
(but above vitrification!) the system remains at a delicate 
balance between energetic and entropic contributions to 
its free energy (see also and the discussion of the 
surface energy below). 

As already mentioned, the RFOT theory has built 
a constructive, microscopic picture of structural relax- 
ations in deeply supercoooled liquids and frozen glasses. 
We learn from the RFOT theory that one may think 
of structural rearrangmeents in deeply supercooled liq- 
uids as activated growth of distinct aperiodic phases 



within each other. The corresponding activation profile 
is UlEI: 



F(N)\ T> T=lVN-Ts c N, 



(22) 



where N is the number of particles in the nucleus of 
the new structural state within the previous structural 
state. The configurational entropy (the ~Ts c N term) 
drives the transitions, while the barrier arises due to the 
mismatch energy penalty, j^/N, between distinct states. 
This mismatch may be thought of as the surface ten- 
sion of the domain wall separating the two alternative 
structural states. As already mentioned, the elemental 
displacements occur on a conditional basis, if an alterna- 
tive structural state is present locally. As a result, the 
full RFOT structural relaxation rate is: 



r RFOT 



2 (T, z = oo)e 



(23) 



(The time scale r str , introduced earlier, was an analog 
of the time scale trfot but in the context of structural 
relaxations in premelted crystalline layers.) 

The elemental translations, corresponding to the 
length d^, play an important role in the RFOT theory, 
for a number of reasons. First of all, the corresponding 
length scale arises as an order parameter during a first 
order cross-over from the "regular" liquid state to the 
one where metastable structures persist for a discernible 
time. (For a deep enough quench, this time is given by 
t rfot f rom Eq.JSBJ).) In the course of this cross-over, the 
localization length of a molecule jumps from infinity to 
a finite value, i.e. d^ itself [2j|. This localization length 
depends only weakl y o n the temperature/density, as the 
DFT study in Ref.|23j suggested. The near constancy 
of the Lindemann ratio d^/a turned out instrumental 
in establishing the near universality of the surface ten- 
sion coefficient between locally competing liquid struc- 
tures in deeply supercooled liquids. This surface tension 
was computed by Xia and Wolynes (XW) without 
adjustable parameters: it depends only logarithmically 
on the Lindemann ratio leading to the following expres- 
sion for the tension coefficient 7: 



2V3vr 



7 



{a/d L f 



Tie 



(24) 



(The notations are from Refs. pH l5q.1 This result, 
among other things, enabled XW to calculate the nu- 
merical value of the barrier for structural reconfiguration 
from Eq. (1221) leading to specific estimates of the size of a 
cooperatively rearranging region. This size grows from a 
few molecular units ("beads"), near T m , to roughly 200 
beads near the glass transition, implying each region is 
about 5-6 beads across near T g . (A bead typically con- 
sists of a few atoms; for a detailed discussion see |57| . 
and also HU). 

The present arguments enable one to compute explic- 
itly the temperature dependence of the Lindemann ratio, 
based on the relation in Ea. i|20[l . The result is shown, 
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again, in FigEU where ( = k B T/{K(T)a 3 (T)) now, see 
Eci. (|21[l . The lattice spacing of the liquid will decrease, 
and the elastic constants will increase, with lowering the 
temperature, albeit weakly. Here of course, we mean 
the high frequency elastic constants, which are definable 
on time scales shorter than the life-times of the long- 
living metastable structures. For most substances, the 
ratio T g /T m ~ 2/3, empirically. As a result, the pa- 
rameter C -1 , in FigO will decrease at most to the value 
of 0.67, before the liquid freezes into a glass, leading to 
d,L(T g ) I 'a(T g ) ~ 0.9, according to the figure. We have 
thus established, on a first principles basis, that the d^j 'a 
ratio indeed varies little with temperature, at most by 
10%. This corroborates the use of the Lindemann crite- 
rion, by the RFOT theory, to establish the near univer- 
sality of the surface tension on the basis of a near univer- 
sality of the molecular displacement at the mechanical 
stability edge, relative to the molecular spacing. Finally, 
the computed temperature dependence of the Lindemann 
length cIl should be measurable by neutron scattering, as 
in Ref.|2g. 

Further we can see explicitly what parameters in the 
problem drive the surface energy 7, by substituting the 
running value of the d^/a ratio from Ea. (|20|l into the 
logarithm in Eq.J33J: 



/csTln 



{a/d L f 



fc B T[ln(4V6) - 1] + V*(T), (25) 



where V t {T) = K{T)a{T)d 2 L {T) is the actual energy bar- 
rier for a typical elemental translation, as stemming from 
the lattice strain; also, [ln(4\/6) — 1] — 1.28. It follows 
from the equation above that the energetic and entropic 
contributions to the surface tension are actually compa- 
rable, at all temperatures above vitrification: T > T g . 
For the reference, we show in Fig0]the V^{T)/kBT ra- 
tio, which is easily computed from Ea. (|20|l . as a function 
of (. At - -67, V t ~ 1.2lk B T. It is not entirely 
clear, at present, what significance should be attributed 
to the £~ ^interval that corresponds to temperatures be- 
low vitrification. 



Finally, it seems instructive to make an explicit es- 
timate of the T~J. apc (T, z = 00). As already mentioned, 
the number of collisions per unit time, for a harmonically 
confined particle, scales linearly with energy, leading to 
(T) / T auto (T m ) = T m /T. This and Eqs.l5|8J yield: 



/378^e 



nddrs 



V*/k B T = T -l (Tm)e l-V*,k BTi 

where D(T m ) and to(T to ) are the regular diffusion con- 
stant and the equilibration time scale tq, of the corre- 
sponding liquid at the melting temperature T m . To(T m ) is 
about ten picoseconds, and the exponential part actually 
does not lead to an activated temperature dependence, 
but a much weaker one, in view of Fig0| 




FIG. 4: Shown is the ratio of the energy barrier to the tem- 
perature, for an elemental translation, as a function of f. 



IV. SUMMARY 

We have established the existence of a universal crite- 
rion of melting, in terms of the ratio of length scales 
characterizing the escape of a particle from its cur- 
rent metastable configuration. The criterion is therefore 
purely kinematic. The obtained quantitative results are 
consistent with earlier studies of displacements at crystal 
surfaces, for several specific substances 0,0,^3- The 
said length scales are closely related to the vibrational 
amplitudes in the crystal bulk, which was argued to un- 
derlie the otherwise puzzling consistency of the empirical 
Lindemann criterion. 

A proper treatment of bulk mechanical stability has 
been performed, and has required considering "fiduciary" 
alternative structural states in the lattice. Such alterna- 
tive states are known to exist in supercooled liquids and 
glasses, and were argued here to exist in a premelted layer 
at the liquid-crystal interface. A proper criterion of me- 
chanical stability was formulated; it stipulates that the 
vibrational molecular displacements be less than the ele- 
mental displacements that would occur during the multi- 
particle structural transitions. We have seen that vanish- 
ing of the shear modulus is not necessary for the lattice 
to be mechanically unstable, consistent with the appar- 
ent high configurational entropy of premelt layers. Direct 
observation of cooperative rearrangements in such layers 
is difficult, but may be possible with available surface 
techniques. 

We have computed the dependence of the elemental 
displacements on the material's stiffness, and on the 
temperature. When taking place in supercooled liquids, 
these can be measured, for instance, by neutron scatter- 
ing. 
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